Modes of magnetic resonance of S=l dimer chain compound NTENP. 
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The spin dynamics of a quasi one dimensional S — 1 bond alternating spin-gap antiferromagnet 
Ni(CgH24N4)N02(ClC>4) (abbreviated as NTENP) is studied by means of electron spin resonance 
(ESR) technique. Five modes of ESR transitions are observed and identified: transitions between 
singlet ground state and excited triplet states, three modes of transitions between spin sublevels of 
collective triplet states and antiferromagnetic resonance absorption in the field-induced antiferro- 
magnetically ordered phase. Singlet-triplet and intra-triplet modes demonstrate a doublet structure 
which is due to two maxima in the density of magnon states in the low-frequency range. A joint 
analysis of the observed spectra and other experimental results allows to test the applicability of 
the fermionic and bosonic models. We conclude that the fermionic approach is more appropriate 
for the particular case of NTENP. 

PACS numbers: 76.30-v, 75.10.Pq, 75.30.Kz 
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INTRODUCTION. 



Dielectric magnetic crystals with an antiferromagnetic 
exchange interaction usually demonstrate Neel ordering 
at some critical temperature. Nevertheless, a number 
of antiferromagnetic systems do not order at temper- 
atures far below the Curie- Weiss temperature. Exam- 
ples are to be found among quasi one dimensional mag- 
netic structures, networks of coupled magnetic dimers, 
frustrated magnets etc. Some of these systems have a 
gap of exchange origin in the magnetic excitation spec- 
trum. This gap stabilizes the spin-disordered ground 
state against perturbations, such as interactions between 
the spin chains, anisotropy, etc. Spin-gap magnets re- 
main in a disordered ("spin- liquid") state down to zero 
temperature. The energy gap in the spectrum of triplet 
excitations may be varied by applying magnetic field or 
pressure. This provides an opportunity to control the 
stability of the spin-liquid state. When the gap vanishes 
at some critical field, the spins order antiferromagncti- 
cally. This phenomenon has been observed and studied 
recently for different spin-gap systems, such as the dimer 
magnet TlCuCl 3 i2, the Haldane magnet PbNi 2 V 2 5 ^, 
the spin 5=1 anisotropic chain system NiCl2-4SC(NH2)2 
(abbreviated as DTN)££. Field induced antiferromag- 
netic ordering in spin-liquid systems was intensively dis- 
cussed, in particular, as a kind of Bosc-Einstcin conden- 



sation of magnons (see, e.g. Rcfs. 

Elementary excitations of a Heisenberg spin-gap mag- 
net are triplets, which may be classified by their spin 
projection S z and momentum k. At zero field the low- 
est excited state is three-fold degenerate (S z = ±1;0). 
The energy of the spin sublevels S z = ±1 depend on 
the magnetic field: E(S Z = ±1) = A =p 9^bH. The 
energy gap vanishes at a critical field H c = A/gfiB- Rel- 
ativistic terms (e.g., crystal field, dipole-dipole interac- 
tion or anisotropic exchange) violate the high symme- 
try of the Heisenberg Hamiltonian. This results in a 
mixing of different spin states and in the lifting of the 
zero-field degeneracy of triplet spin sublevels. In this 
case the field dependence of spin sublevels, and hence 
the relation between the zero field energy and the criti- 
cal field, become more complicated. There are three ap- 
proaches to describing the spectrum in a magnetic field 
and calculating the critical field: fermionic models*^, 
the bosonic (or macroscopic) approac h 10 ' 11 , and pertur- 
bative calculation s 12 ' 13 . The fermionic approach is jus- 
tified for the case of a purely one-dimensional spin-gap 
system. The bosonic approach implies proximity of three 
dimensional magnetic order. It is valid if the energy gap 
is much smaller then the main exchange integral. The 
perturbative approach is valid at low fields g^sH <C A 
irrespectively to the system dimensionality. It demon- 
strates that quantized collective excitations are equiva- 
lent to sublevels of a single spin S — 1 in an effective 
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crystal field. Predictions of the perturbative model co- 
incide with those of fermionic model (see Ref.12;), thus 
expanding applicability of the perturbative model in the 
one-dimensional case up to critical field. If the zero field 
spin sublevels are fixed, both types of results are indis- 
tinguishable at low field. Nevertheless, they differ in the 
vicinity of critical field: the perturbative/fermionic ap- 
proach predicts a closing of the gap linearly with the mag- 
netic field, while macroscopic approach predicts a faster 
square root dependence. The scenario of the spin gap be- 
haviour in a magnetic field is not universal: the Haldanc 
compound NENP was found to fit better to the fermionic 
models, while another Haldane magnet PbNi2V2C>8 and 
the dimer-type spin-liquid magnet TICUCI3 fit to the 
macroscopic model^ii. Thus, it is of importance to study 
further spin-gap systems which demonstrate closing of a 
spin gap in a magnetic field. 

Another problem that we address in this paper is 
the stability of excitations in spin liquid systems. The 
magnetic excitations spectrum usually includes several 
branches which may form two-particle continua. A de- 
cay of magnons is possible, if a single-magnon branch 
enters such a continuum. In this case magnon modes be- 
come overdamped, often resulting in a termination of the 
single-particle spectrum. 

The recently found compound 

Ni(C 9 H24N 4 )N02(C104), abbreviated as NTENP, 
is a good test probe for spin excitations at a critical 
field closing the gap, as well as for field induced ordering 
and spectrum termination^^— Besides, it belongs to a 
new family of spin-gap magnets with dimerised S = 1 
spin chains. In crystals of NTENP, Ni 2+ magnetic 
ions form dimerised spin S=l chains running along a 
axis of the triclinic crystal. In-chain exchange integral 
takes two alternating values J\ and J2 with a ratio 
a = J1/J2 = 0.45. The larger exchange integral is 
h/kB = 54.2 KJ£ The spin gap in S — 1 chains may be 
of a Haldane type, if chains arc uniform (a = I), or of a 
dimer nature in the opposite limit case a = 0. There is 
a zero gap point between these limiting cases at a ~0.6 
(see, e.g., Ref. HI), thus NTENP belongs to the dimer 
class. The ground state of NTENP is a nonmagnetic sin- 
glet separated from the triplet excitations by an energy 
gap of about 1 meV. The triplet energy levels are split 
by the crystal field: the energies of the sublevels, found 
in neutron scattering experiments are 1.06, 1.15 and 1.96 
meV^ respectively. In zero magnetic field, the doubled 
energy of the lowest sublevel is a little bit lower than the 
energy of the upper sublevel. The excitations energies 
arc field-dependent, and in a magnetic field the decay of 
magnons becomes allowed by the energy conservation 
lawi 16 i 18 The critical field H c , closing the spin-gap is 
anisotropic: it varies from 85 kOe in a magnetic field 
applied along the chains to 114 kOe in the perpendicular 
orientation^ The field-induced ordering takes place 
only at low temperatures: in a magnetic field of 100 kOe, 
applied along the chains, the Neel temperature is about 
0.5 K. The phase diagrams for the induced ordering at 



different orientations of the magnetic field are described 
in Refs. [U [22|. Antiferromagnetic ordering above H c is 
confirmed by means of neutron scattering^ 

In the present paper we report the detailed electron 
spin resonance (ESR) study of this material. This tech- 
nique allowed us to detect excitations with energies 
of about 0.05meV (corresponding to the resonance fre- 
quency of approximately 10 GHz) with the resolution of 
about 0.005mcV (1 GHz). We observed magnetic res- 
onance modes of thermally activated triplet excitations 
and the over-gap transitions, as well as antiferromagnetic 
resonance absorption in the field-induced ordered phase. 
The splitting of triplet modes was observed, which can be 
ascribed to the dispersion of the triplet excitations with 
wave vectors along the direction of a weak exchange (i.e. 
perpendicular to chains). Further, the spin resonance 
data provided a possibility for comparison of fermionic 
(perturbative) and bosonic (macroscopic) models. We 
found that our data favour the fermionic description. We 
refined parameters of the effective Hamiltonian of the 
triplet excitations and estimate the inter-chain exchange 
integral. The effects, related two the two-magnons decay 
of spin excitations, were observed. 



II. EXPERIMENTAL DETAILS AND SAMPLES. 

A set of custom-built transmission type ESR spectrom- 
eters covering the frequency range from 9 to 260 GHz 
was used for the experiments. The spectrometers were 
inserted in a 4 He-bath cryostats equipped with 60 — 120 
kOe cryomagnets. Microwave radiation was produced by 
commercial Gunn diode generators and backward wave 
oscillators. The microwave resonator with the sample in- 
side was placed in a hermetically sealed chamber filled 
with heat exchange gas. The temperature of the sample 
may be stabilized by means of a heater in the interval 
between 1.3 K and 30 K. Experiments in the tempera- 
ture range 0.4 — 1 K were performed in a custom-built 
spectrometer equipped with a 3 He-3 pumping refrigera- 
tor. Microwave modes of several resonators were used to 
cover the entire frequency range. 

In an experiment the signal transmitted through the 
resonator is recorded as a function of magnetic field. A 
diminishing of transmission means microwave absorption. 
Resonance absorption appears if the microwave frequency 
equals the frequency of a transition between states \i), 
\j) {fiuj — \E tJ \ = \E t — Ej\). The following selection 
rules should be satisfied: (i) the dipolar transition matrix 
element («|(hS)|j) is non-zero (here h is the amplitude of 
the oscillating field) and (ii) states \i) and \j) have the 
same wavevector k. The total ESR absorption is a sum 
over all pairs of spin states and an integral over the k- 
space. 

The configuration of the microwave field in the res- 
onator is different for different microwave modes, hence 
the polarization of microwave magnetic field in the sam- 
ple is also different. Due to the finite sample size, even for 
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FIG. 1. (color online) ESR absorption above IK at H||c* (left) 
and H||o (right). Vertical dashed lines mark the position of 
the paramagnetic resonance with g — 2.00. 



a given mode, different polarizations of microwave field 
are present within the sample. Therefore, for the typical 
sample size, and for most of microwave resonator modes, 
the polarization selection rules do not prohibit any modes 
of magnetic resonance absorption. 

We have used the samples of deuterated NTENP from 
the same batch as in RefJl7l. 
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FIG. 2. (color online) Temperature dependences of the in- 
tegrated intensity of ESR absorption. The dashed straight 
line on the main panel corresponds to the exponential law 
I = Aexp{ — ^} with E = 10. 4K. Inset: Temperature de- 
pendence of half-linewidth. The dashed line at the inset is a 
guide to the eye. 



III. EXPERIMENTAL RESULTS. 

A. Above IK: ESR of thermally activated 
excitations. 

The intensity of the magnetic resonance absorption ob- 
served above IK increases with temperature up to about 
20 K. Upon heating beyond 20 K the lines broaden and 
disappear. Examples of the temperature evolution of spin 
resonance absorption for two orientations of the mag- 
netic field are presented in Fig. [TJ The lines are strongly 
shifted from the paramagnetic resonance position of free 
electron spins. The resonance field is anisotropic, which 
is in agreement with the zero-field splitting of the triplet 
sublcvels observed in neutron scattering experiments. 
Below 5K the integrated intensity follows a law of thermal 
activation / ~ exp{— Eq/T} with the activation energy 
Eq = 10.4 ± 1.0K (see Figure [5]). This activation energy 
corresponds well to the known energy of lowest excited 
state of 1.06 meVJ£ Exponential freezing out of inten- 
sity proves that this signal is due to transitions between 
sublcvels of excited states. The temperature dependence 
of the linewidth is shown on the inset of Figure [2]) . The 
linewidth at T=4K is practically the same as in the low- 
temperature limit. The resonance field at 4K also coin- 
cides with its low-temperature value. Thus, absorption 
at T=4K may be attributed to noninteracting excitations 
distributed near the bottom of the spectrum. 



ESR absorption lines were collected at a set of fre- 
quencies at T=4.2 K for three mutually orthogonal ori- 
entations: H||a, H||c* and H _L a,c*. Here a is the 
crystallographic axis of the triclinic crystal and c* is per- 
pendicular to the crystallographic (afe)-planc. From here 
on we will mark the resonance modes observed in the 
high temperature range T >1 K as follows: HT1, HT2 
and HT3 are the modes observed at H||a, numbered in 
the order of ascend of zero-field frequency , HT4 is the 
mode observed at H||c* and HT5 is the mode observed 
at H _L a, c* . Examples of the resonance lines at different 
frequencies are given in Figure [3] 

At H| \a we observe single resonance lines for all modes. 
For H||c* or H _L a, c* at low temperature the absorp- 
tion line is split into two close components (see Figure 
[21). The intensities of the split components are approxi- 
mately equal. The high frequency modes HT2 and HT3, 
observed in the H||a orientation, are much broader then 
other resonance modes. 



B. Below IK: Singlet-triplet transition and 
antiferromagnetic resonance. 

The experiments at T < 1 K were performed only for 
H||a because the critical field H c is a minimum at this 
orientation, enabling the largest value of H/H c . 

At temperatures below 1 K the triplet excitations are 
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FIG. 3. (color online) ESR absorption lines for different 
resonance modes. All absorption curves, unless specified oth- 
erwise, are measured at 4.2K. HT1...HT5 denote resonance 
modes as described in the text. The unlabelled narrow ab- 
sorption line is a DPPH marker (g — 2.00). 



frozen out and do not contribute to resonance absorption. 
Nevertheless, other types of the ESR absorption arise at 
high fields (Figure 0]). First, two absorption components 
(LT1 and LT2) split by about 10 kOe are observed be- 
low 90 kOe. These resonances shift to lower fields with 
the frequency increase. The second low-temperature ab- 
sorption range is found above 90 kOe. Here two broad 
bands of microwave absorption denoted as LT3 and LT4 
are found. Both of the absorption maxima shift to the 
higher fields with increasing frequency, as shown on Fig- 
ure [5] 

Mode LT1 is usually much weaker then LT2. The rel- 
ative intensity of these modes varies irregularly with the 
frequency. This behaviour is, probably, due to the depen- 
dence of their excitation on the dominating polarization 
of microwave field. In contrast to HT-modes, intensities 
of modes LT1 and LT2 decrease with increasing tempera- 
ture (see Figure O, i.e. they are not thermally activated. 
These frequency-field and temperature dependences al- 
low us to conclude that modes LT1 and LT2 correspond 
to singlet-triplet transitions. Two other low-temperature 
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FIG. 4. (color online) ESR absorption at 0.4K at different fre- 
quencies. LT1, LT2, LT3, LT4 are the absorption components 
(see text).H||a 
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FIG. 5. (color online) Frequency-field diagram for modes ob- 
served at T — 0.4K. Lines are guides to the eye. H||a 



modes (LT3 and LT4) correspond to resonance modes of 
the field-induced ordered phase (above H c ) and may be 
considered as antifcrromagnctic resonance modes. 

Upon heating above IK, the resonances LT1 and LT2 
are no longer resolved. At several frequencies around 
80 — 100 GHz we observed an indication of weak ESR 
mode at a temperature of 1.2K. This mode is, probably, 
related to modes LT1 and LT2. It is not split and located 
between LT1 and LT2 modes. A weak trace of the critical 
region survives to about 2.5 K in the form of a wide res- 
onance, see Figure [5J We do not observe a clear anomaly 
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FIG. 6. (color online) Temperature evolution of low- 
temperature modes at / =35.25 GHz, H || a. The vertical 
dashed line marks the position of a paramagnetic resonance 
with #=2.00. 
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FIG. 7. (color online) Temperature evolution of low- 
temperature modes through the Neel temperature at 
/ =28.25 GHz, H||a. Gray curves — experiment, thick curves 
(red) — guides to the eye obtained as a two-lorentzian best 
fit. 



in the temperature evolution of LT3 and LT4 modes at 
the Neel temperature (which is 0.6 K at the field of 100 
kO o 21 i 22 ) . Instead, these modes merge with the wide ab- 
sorption near the critical field. This transformation is 
shown on Figure [7] 
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FIG. 8. (color online) Frequency-field diagram, showing all 
observed resonance modes (symbols). Curves show different 
model approaches discussed in the text. Dotted line : macro- 
scopic (bosonic) model. Solid and dashed lines: perturbative 
(fermionic) model in different points of k-space. 



C. Frequency-field diagrams. 



Resonance fields measured in different orientations are 
collected on the frequency-field diagram (Figure [8]). We 
observe three thermally activated modes in H||a orienta- 
tion (HT1, HT2, HT3), one thermally activated mode 
at H||c* (HT4) and one thermally activated mode at 
H _L a,c* (HT5). Positions of these modes are shown 
as measured at 4.2 K. At T =0.4 K, in the H||a orienta- 
tion, we observe four modes (LT1, LT2, LT3, LT4). The 
frequencies of these modes are shown as measured at 0.4 
K. A weak mode measured at 1.2 K originates, probably, 
from the temperature evolution of modes LT1, LT2 and 
is also plotted here. The splitting of modes HT4 and 
HT5 described above is too small to be resolvable on the 
Figure |H1 being within the error bars. 

Thermally activated modes HT1, HT2 and HT3 
demonstrate zero-field gaps of 20 ± 3, 190 ± 2 and 
210 ± 2 GHz correspondingly. The frequencies of LT1 
and LT2 modes may be approximated with a linear law 
/ = -7/J+A with 7 = 2.8GHz/kOe (assuming g = 2.00) 
and a gap A of 235 ± 5 and 265 ± 8 GHz for modes LT1 
and LT2 respectively. Critical fields for these modes, 
estimated as H c = A/7, are 84 ± 2 and 94 ± 3 kOe re- 
spectively. 
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FIG. 9. (color online) Neutron scattering data on the field 
dependence of triplet sublevels' energies from Refs[l^, llfil . 
Curves: predictions of the perturbative model (solid) and 
macroscopic (dashed) model. H||6 



IV. DISCUSSION. 

This section is organized as follows: First, we will 
briefly summarize our experimental data and compare 
them with earlier results. Second, we will check the cor- 
respondence of perturbative and macroscopic approaches 
of spectrum description to experimental data. Finally, we 
will estimate the inter-chain coupling from the observed 
splitting of resonance modes and will analyze possible in- 
fluence of magnon decay on the spectrum terminating in 
NTENP. 



A. Qualitative analysis. 

Three types of spin resonance modes were observed: 
thermally activated intra-triplct transitions, singlet- 
triplet transitions and antiferromagnetic resonance ab- 
sorption. 

The zero field singlet-triplet transition frequencies es- 
timated from the linear approximation of modes LT1 and 
LT2 (235 ± 5 and 265 ± 8 GHz, correspondingly) are in 
a good agreement with the lowest triplet excitation ob- 
served in neutron scattering experiments^ at 235±5GHz 
(i.e. 0.97±0.02mcV). The frequencies of two other excita- 
tions were determined from neutron scattering as 273 ±5 
and 484 ± 24 GHz (1.13 ± 0.02 and 2.0 ± 0.1 meV). The 
zero-field ESR intratriplet transitions at 20 ± 3, 190 ± 2 
and 210 ± 2 GHz correspond well to the triplet sublevels 
energy differences, measured in neutron scattering exper- 
iments as 17 ± 7, 211 ± 25 and 249 ± 25 GHz. 



The derivative -^j for the HT1 mode is approximately 
twice as large as expected for a free spins resonance with 
a g-factor of 2.00. Besides, the zero field frequency of 
the HT1 mode is much smaller than the splitting of the 
triplet levels due to the main anisotropy term. This fea- 
ture is well known for spin S = 1 as a " two- quantum" 
transition between the states with S z = +1 and S z = — 1 
(see, e.g. Refl23h. This type of the /(-ff)-dependence 
should be observed for the direction of magnetic field 
along the main anisotropy axis and is very sensitive to 
the field orientation. Thus, our observation of HT1 mode 
proves that a is close to the main anisotropy axis. 

The extrapolation of modes LT1 and LT2 to zero fre- 
quency results in the critical fields of 84 ± 2 and 94 ± 3 
kOe. This is again in agreement with the magnetization 
measurements^ on non-deuterated sample of NTENP 
done at 1.3 K where the increase of magnetization starts 
at 84±1 kOe for H j| a. A higher value of the critical 
field H c (H\\a) = 91.7±0.1kOe was reported^ for deuter- 
ated NTENP. This value was extracted from the critical 
exponent fit of the phase boundary of the field-induced 
ordering. However, due to the lack of the experimental 
data points on the phase diagram in the immediate vicin- 
ity of the critical field this estimate of the H c should be 
considered with a higher error bar. 

The singlet-triplet transitions (mode LT1 and LT2) 
should be forbidden in the Heiscnbcrg case as transitions 
between different spin multiplets. " Two- quantum" tran- 
sitions (mode HT1) are also forbidden in the case of axial 
anisotropy, because there is no dipolar matrix element 
connecting states with S z = 1 and S z = —1. In the 
presence of in-plane anisotropy it becomes allowed and 
should be observable at microwave magnetic field polar- 
ized parallel to static field. The modes HT2 and HT3 
correspond to allowed transitions with large matrix ele- 
ments at a conventional ESR polarization ( AS Z = il)- 
Despite this, the modes HT2,3 have weaker intensities 
and are broader than "forbidden" modes. The intensity 
and lineshape for different modes will be analyzed in Sub- 
section D. 

Finally, there is a splitting of certain resonance modes: 
(i) Thermally activated modes HT4 and HT5 demon- 
strate splitting into two components of approximately 
equal intensity, (ii) Instead of a single-mode singlet- 
triplet transition, we observe two modes LT1 and LT2, 
LT1 is usually much weaker. This strong difference of 
LT1,2 modes intensities was repeated regularly for sev- 
eral samples. This excludes possible explanation of the 
splitting being caused by crystallites of different orienta- 
tions. 



B. Comparison of the results of perturbative 
(fermionic) and macroscopic (bosonic) approaches to 
the low-energy dynamics of NTENP 

We fit the observed ESR frequencies to the pertur- 
bative and macroscopic models mentioned in the Intro- 
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TABLE I. Comparison of perturbative (fermionic) and macro- 
scopic (bosonic) models with the known experimental results. 
Parameters of the models are obtained from the fits to the 
ESR data presented in Figure [8] 
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includes 5 parameters: A, D, E, g and an angle <\> be- 
tween anisotropy x axis and c* direction. Given these 
parameters, eigenenergies of the Hamiltonian (JTJ) can be 
found in any orientation of magnetic field. Differences 
between the eigenenergies correspond to frequencies of 
thermally activated modes. The lowest eigenenergy cor- 
responds to the singlet-triplet transition. The zero-field 
energies for this model are A, A + D + E and A + D — E. 
The field dependences can be easily found for H||z: 

E ia = A + D±^(gii B H) 2 + £ 2 (2) 
E 3 = A (3) 

For an arbitrary direction an analytical solution be- 
comes too cumbersome. Critical fields for H||z and 
Hlz are 



duction. The k-dependence of model parameters will be 
neglected in this Subsection because at low temperatures 
the quasiparticles are excited at the bottom of the spec- 
trum. Consideration of the k-depcndcncc of energy levels 
near the bottom of the spectrum results in the asymmet- 
ric distortion of the absorption line, while the maximum 
of absorption still corresponds to the bottom of the spec- 
trum due to a large statistical weight of the corresponding 
states^. The main anisotropy axis z is taken to be par- 
allel to the chain direction a, this choice is based on the 
observation of the "two-quantum" transition mode HT1. 
The second anisotropy axis may not be chosen from the 
symmetry consideration of the triclinic crystal and should 
be derived from the experiment. Thus, the direction of 
the second anisotropy axes in a plane perpendicular to a 
is fitting parameter of the model. Finally, we will assume 
^-factor to be isotropic. 

The perturbative mode l 12 ' 13 treats triplet excitations 
as non- interacting S—l quasiparticles in the effective 
crystal field and in the external magnetic field. This ap- 
proach is limited to low fields when the Zeeman energy 
is small in comparison with the spin gap. However, as it 
was shown in Reffl2l. in a one dimensional case the results 
of perturbative formalism coincides with the fermionic 
formalism of ReflU, which is valid at all fields. For this 
reason, taking into account the strongly one-dimensional 
character of spin interactions in NTENP, we will apply 
the results of the perturbative approach all the way up 
to the critical field. 

The effective Hamiltonian of the perturbative model 

is: 



H 



tripl 



D S 



E S 



S y 



-g^B (HS) 



(i) 



here A is the gap separating singlet state with zero en- 
ergy from the excited triplets, D and E are the effective 
anisotropy constants, and g is the ^-factor. This model 



(gHBH c 



(A + D) - E 



2 

1^2 



E 2 



A + D — E cos 2a 



(4) 
(5) 



here a is the angle in (xy) plane counted from x direction. 

The model was fit to the experimental data 
for thermally activated intertriplet transitions (modes 
HT1, 2, 3, 4, 5) and for LT1 as a singlet-triplet tran- 
sition. The best fit corresponds to the following model 
parameters: D = -197 ± 5 GHz, E = 10.5 ± 1.0 
GHz, A = 434 ± 5 GHz, g = 2.02 ± 0.05 and <$> = 
(0.59 ± 0.04)rad= (34 ± 2)°. The zero-field energies of 
triplet sublevels in frequency units, as predicted by this 
model, are collected in Table HI The modelled frequency- 
field dependences are presented in Figure [8] by solid lines. 
They give a good approximation for the observed reso- 
nance modes. 

The macroscopic (or bosonic) approac h 10 ' 11 describes 
the low energy spin dynamic of the spin-gap magnet as 
oscillations of a vector field 77. This model assumes that 
the spin oscillation frequencies are small with respect to 
the exchange frequency and that antifcrromagnetic or- 
dering appears above the critical field. Thus, the model 
requires non-zero inter-chain coupling. 

The Lagrangian of the macroscopic model (without 
gradient terms) is 



L 



1 A 
-(77 + 7 [r/xH]) 2 --7, 2 + 

+ b j(vl + V 2 y -2vl)+ b -§(r) 2 x -V 2 y ). 



(6) 



here A is the exchange constant describing the energy 
gap, bi and hi are the effective anisotropy constants, 
7 is a gyromagnetic ratio (7 = (,9/2.0)70, where 70 = 
2.80GHz/kOe is a free electron gyromagnetic ratio). Be- 
low H c , the corresponding dynamic equation is: 
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The effect of inter-chain interaction. 





-7 2 H(r? 


■ H) + At] — 




( V, 




( V, \ 


bi 


Vy 


j -b 2 


-Vy = 0. (7) 




{-2V* 




^ / 



For the case of H||z the excitation energies are: 



E l,2 



r-H 1 +A — bi± 



±J4- f 2 H 2 {A-b 1 )+bl 



E\ = A + 2h 



The critical fields for Hllz and HI 2 arc 



7 2 iJ 2 N = A-h 



7 2 # c 2 1 = mm{A + 2b 1 ; 



A — b\ — 62 cos 2a 



(8) 
(9) 



(10) 
(11) 



Note, that for b\ > (i.e. Eyp < E3, which corre- 
sponds to the case of NTENP) this model predicts no 
anisotropy of critical field for 62 = 0. 

This model also has 5 fit parameters: A, b\, 62, 7 
and the orientation of the x axis with respect to the 
c* direction </>. The best fit obtained is shown on 
the Figure [H It corresponds to the following parame- 
ters: A = 96700 ± 2000GHz 2 , b x = 43600 ± 1000GHz 2 , 
b 2 = 4640 ± 500GHz 2 , 7 = 2.71 ± O.lOGHz/kOe (i.e. 
g = 1.94), (f> = (0.61 ± 0.08)rad= (35 ± 4)°. The critical 
fields and zero-field energies of triplet states, predicted 
by this model, are collected in Table fl] This macroscopic 
model also yields predictions close to the experimental 
data, but the agreement is less satisfactory. In particu- 
lar, systematic deviations arc found for modes HT1, 4, 5. 

We conclude that the ESR data (Figure [8]) favours the 
perturbative model though not unambiguously. To pro- 
vide an additional check, we have calculated field depen- 
dences of energy levels in the H||6 orientation of mag- 
netic field using the model parameters values determined 
above. These model predictions can be compared with 
independent data from neutron scattering, obtained in 
the same orientation of magnetic fiel d 16 ' 18 . As born out 
in Figure GO predictions of the macroscopic model visibly 
deviate from the experimental points, while the pertur- 
bative model describe the neutron scattering data fairly 
well. 

Thus we conclude that the perturbative (fcrmionic) 
model gives a better correspondence with the exper- 
iments for the case of NTENP. Disagreement of the 
macroscopic model with the experiment is probably due 
to the violation of low-frequency condition, since the en- 
ergy of the upper triplet component of 450 GHz is about 
40% of the largest exchange J%/h. 



Weak inter-chain interactions estimated as J± ~ 
0.01meV« 0.1K were reported in Refli~7l. These interac- 
tions will result in a dispersion of excitations propagating 
transverse to the chains. To take into account inter-chain 
coupling in the model calculations we will use form of 
excitations spectrum that combines the energy gaps and 
dispersive terms analogously to an antiferromagnet with 
anisotropic spin wave velocities: 



E = A 



-2s 



|S_L 



2 sin 2 (ka) + s\ sin 2 (kb) 
(1 — cos(ka) cos(kb)) 



+ 



(12) 



here s\\ t ±_ = 2SJ\\^±_ are spin-wave stiffnesses along the 
chain and perpendicular to the chain, correspondingly 
(we use exchange energy in the form E = j\ JijSiSj), 
a and b arc translation vectors of the rectangular lattice 
along the chains and in the transverse direction respec- 
tively, k is the vector in k-space relative to the point of 
minimum of energy, and A; (i = 1,2,3) are the corre- 
sponding gaps, respectively. 

The longitudinal spin- waves stiffness si is due to the 
dominating in-chain exchange interaction and has a value 
o£il: S|| = 8.6meV= 2080GHz. For a weak inter-chain 
interaction s± <C sm, and the dispersion for excitations 
propagating along the 6-axis is: 



E-, 



ka=0 



2S||S_l(1 



"(I 



cos(kb)) 
cos(kb)) 



(13) 



From the exchange constants obtained in neutron scat- 
tering experiments we have j} B ^. ~ 1.6K. Thus, all exci- 
tations with transverse wavevectors should be populated 
above IK, including those at kb = n. Therefore, both 
points with maximum density of states of excitations will 
contribute to microwave absorption: the bottom of the 
spectrum at kb = and the saddle point at kb = ir. The 
parameters of effective Hamiltonian (TTJ) in these points 
(Ao, Do, Eq and A n , D T , E w , respectively) are slightly 
different. This will result in the splitting of the resonance 
modes. 

The two sets of effective Hamiltonian parameters are 
related as follows (see Eqns. ©, ©, (fT3")l ): 



A, 



Aq + 4s||Sj_ 



v /(Ao + £>o + g )^+4s||sI 
2 

V /(A + A)-£o) 2 + 4s||/ 
2 

v /(A + £> + ^ ) 2 +4s||s7 
2 

v /(A +D - J E ) 2 + 4s||S_ 



(14) 



(15) 



A. (16) 
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Now we can ascribe LT1 and LT2 modes to the transi- 
tions between the ground state and kb = and kb = 7r 
points in the excitations spectrum. This allows us to 
extract the value of (s||Sj_) from the splitting of these 
modes. The best fit corresponds to ^/s^s± = 60 ± 5GHz 
(0.25meV). 

The inter-chain exchange integral can be estimated as 
follows: Ji = (sj_/s||)J|| = ((sxS||)/sfj)J|| - 0.05 K. 
This value is quite small and there are other interac- 
tions (i.e. dipolar, or anisotropic exchange interactions) 
of the same strength. Thus, the appearance of a weak 
transverse dispersion may be due to the combined action 
of several interactions, which are much weaker than the 
dominating in-chain exchange. 

Frequency-field dependences for both the kb = and 
kb = 7r excitations, calculated for all modes, are shown 
for comparison on the Figure \8\ The splitting of reso- 
nance lines are well described by our dispersive model: 
The splitting is largest for the singlet-triplet transition 
(modes LT1 and LT2). It is due to the different values 
of lowest triplet sublevel energy in different points of k- 
space. For the thermally activated modes, the splitting is 
much weaker, since it is determined by the difference of 
the triplet sublevels energies and all sublevels are shifted 
in the same direction due to the dispersion. For the HT1 
mode the splitting is not observable because of the steep 
increase of frequency with magnetic field, while for modes 
HT4 and HT5 a weaker frequency-field dependence al- 
lows us to resolve the doublet by scanning the field. The 
predicted splitting of modes HT4 and HT5 is about IkOe 
which is in agreement with experimental observations. 
As for the high-frequency modes HT2 and HT3 (H||a), 
the predicted splitting of approximately 5kOc may not be 
observed because of the large linewidths of these modes 
(about 5kOc). 

Equal intensities of the components of the split doublet 
for thermally activated modes HT4 and HT5 are due to 
almost equal population numbers at the top and bottom 
of the spectrum of transversal excitations. The obser- 
vation of both LT1 and LT2 modes indicates that the 
ground state is a mixture of kb = and kb = n states. 

D. Indications of the quasiparticles decay. 

Taking into account the real values of the derivative 
df/dH, the linewidth of modes HT2, HT3 is about 20 
GHz, while the linewidth of HT1 is 1 GHz. For modes 
HT4,5 it is about 4 GHz (splitting included). 

We performed a modelling of the ESR lineshape by 
integrating the expected signals over the k-space. 

Model calculations were performed in the perturbative 
approach using the best fit parameter values obtained 
as described in previous sections. The parameters Ak, 
Dk and £"k for the effective Hamiltonian (jTJ) where cal- 
culated using the zero field energies given in Table U for 
the perturbative model and using Ean.(|12[). Then eigen- 
frcquencics and cigenfunctions for Hamiltonian (TTJ) and 




20 40 



H [kOe] 

FIG. 10. (colour online) Modelled lineshapes of ESR absorp- 
tion at 4K for H| \a. Solid lines: modelling taking into account 
dispersion both along the chain and in the transverse direc- 
tion. Dashed lines: modelling taking into account only dis- 
persion along the chains. Solid line with symbols at 210GHz: 
modelling taking into account higher inherent linewidth of the 
Sz = magnon. Dashed lines with symbols: guides to the eye 
marking position of different modes at different frequencies. 



transition matrix elements were found. Finally, we cal- 
culated absorption by summing over the k-spacc, tak- 
ing into account the population of the sublevels involved 
in a given transition. The lineshape of individual tran- 
sition was assumed to be Lorentzian with the inherent 
width of 0.5 GHz. This value was found to reproduce 
well the low-temperature linewidth of the narrow HT1 
mode. The inherent linewidth was assumed to be the 
same for all transitions. The absorption was calculated 
for optimal polarization of the microwave magnetic field, 
i.e. for HT1 mode the longitudinal polarization of mi- 
crowave field was used, while HT2 and HT3 modes ab- 
sorption was calculated at transverse polarization. In an 
experiment both polarizations are usually present. For 
low frequency modes, the transverse polarization domi- 
nates. 

Results of this modelling are presented on Figure [TO] 
The model predicts that the absorption by modes HT2, 
HT3 should be in a wide band with a sharp edge corre- 
sponding to kb = and a peak corresponding to kb = 7r, 
instead of the wide and approximately Lorentzian lines 
observed. The model predicts also, that the resonance 
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20 40 60 80 

H [kOe] 



FIG. 11. (colour online) Modelled field dependences of triplet 
energy levels (perturbative model). Boundaries of the zone 
due to the dispersion of excitations in the transverse directions 
are marked by solid (kb = 0) and dashed (kb = n) lines. 
Arrows marks transitions corresponding to the observed ESR 
modes. For modes HT1...HT3 only transitions between q±b — 
states are shown for the sake of simplicity. The grey shaded 
area shows the two-particle continuum, corresponding to the 
lowest excitation. 



absorption at the HT2 mode should be large enough to 
be detected at least at frequencies above 100 GHz. It 
disagrees with experimental observations, because, in the 
frequency range 170-240 GHz, a smeared peak is observed 
experimentally. Besides, below this frequency range, the 
resonance HT2 was not observed, while HT1 is observ- 
able. 

A simple way to reconcile the model with the exper- 
iment is to suppose that the inherent linewidth of HT2 
and HT3 modes is much larger then that for HT1 mode. 
The broad modes HT2 and HT3 involve transitions to 
the highest excitations branch, corresponding to S g = 0, 
while the narrow modes HT1, HT4, HT5, LT1, LT2 cor- 
respond to transitions between lower energy branches. 
Therefore, the existence of a strong relaxation specifi- 
cally for the S z = excitations could explain observed 
discrepancy. We model this strong relaxation by assum- 
ing inherent linewidth of HT2 and HT3 modes to be equal 
to 5 GHz instead of 0.5 GHz. This, indeed, results in the 
smeared lincshape of HT3 resonance, as shown in FigfTUl 
in a solid line connecting filled circles. 

Possible reason for strong relaxation for the S z = 
excitations is the entry of this branch into the two parti- 
cle continuum range. Similar explanation was proposed 
in Reflill to describe the disappearance of the highest- 
energy excitations observed in neutron scattering exper- 



iment. The frequency corresponding to the bottom of 
the two particle continuum may be derived as a doubled 
frequency of the k = excitation (for the excitations 
with a positive dispersion). Energies of excitations and 
two-magnon continuum related to the lowest branch are 
shown in Figure QTJ The decays of highest excitation 
branch into two excitations from the lowest branch may 
occur in NTENP due to the close values of the energy of 
highest excitation mode (which are 435 GHz for kb = 
and 450 GHz for kb = it) and the bottom of two particle 
continuum (which is 453GHz for H = 0). In a magnetic 
field, the energy of the lowest branch decreases. At fields 
above 6 kOe the excitations of the highest branch became 
unstable with respect to the decays into two magnons of 
the lower branch. 

It is an open question whether this intersection with 
the continuum destroys the excitation branch or just 
makes the magnon life time shorter. In any case, decay 
should result in a broadening of ESR line. The observed 
disappearance of the HT2 resonance in fields above 20 
kOe (or in the frequency range below 150 GHz) may be 
due to such decays into two particles of lower branch. 



E. Antiferromagnetic resonance 

The antiferromagnetic resonance absorption above H c 
consists of a broad band with two absorption maxima. 
This observation corresponds to the reported extreme 
broadening of the lowest energy antiferromagnetic mode 
in neutron scattering above H c for H||6ii£ The splitting 
of the antiferromagnetic mode into two close resonances 
LT3 and LT4 may be due to the a dispersion in the trans- 
verse direction, analogous to the splitting of LT1,2 and 
HT4,5 modes described above. 



V. CONCLUSIONS. 

The low-energy spin dynamics of the dimcrised S=l 
spin-chain compound NTENP is studied by means of 
ESR. Modes of magnetic resonance corresponding to 
transition between the spin sublevels of the collective 
triplet excitations of the spin-gap magnet are observed 
and identified. The singlet-triplet transitions to a spin 
gap state are also observed. The data analysis demon- 
strates that low-energy spin dynamics of this spin-gap 
magnet may be described by a perturbative model in the 
whole field range up to the critical field. The singlet- 
triplet transition and the intratriplct transitions have a 
doublet structure due to the dispersion of excitations 
propagating perpendicular to the chains. The energy of 
the inter-chain interaction is determined from the dou- 
blet splitting. 

In addition we have observed a strong broadening and 
termination of spin resonance modes at the boundary 
of the two particle continuum, indicating the influence 
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of the two-magnon decay processes on the spectrum of 
excitations in NTENP. 
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